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Abstract 

We have searched for a CP violation signature arising from an electric dipole mo- 
ment {d-r) of the r lepton in the e'^e" — t'^t~ reaction. Using an optimal ob- 
servable method and 29.5 fb~^ of data collected with the Belle detector at the 
KEKB collider at ^ = 10.58 GeV, we find Re{dr) = (1.15 ± 1.70) x lO^^^ecm 
and Im{dr) = (—0.83 ± 0.86) x 10~^^ecm and set the 95% confidence level limits 
-2.2 < Re{dr) < 4.5 (IQ-^'^ecm) and -2.5 < Im{dr) < 0.8 (lO^^'^ecm). 

Key words: tau, electic dipole moment, CP violation, optimal observable 
PACS: 13.40.Gp, 13.35.Dx, 14.60.Fg 
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1 Introduction 



While large CP violating asymmetry has recently been confirmed in B-meson 
decay [1,2], the Standard Model (SM) does not predict any appreciable CP vio- 
lation (CPV) in the lepton sector. However, physics beyond the SM could pro- 
duce CPV in leptonic processes; we would expect such effects to be enhanced 
for r leptons due to their large mass. Parameterizing CPV in r-pair production 
by an electric dipole moment dr, several authors have found possible effects 
of order {d^-l ~ lO^^^ecm due to new physics models [3,4]. The best existing 
bounds on dr are indirect, requiring \dr\ < O(10~^'')ecm based on naturalness 
arguments [5] and precision LEP data [6] .0 Previous direct measurements have 
been performed at LEP, where L3 [9] found —3.1 < Re{dr) < 3.1 x 10~^^ecm 
and OPAL [10] found \Re{dr)\ < 3.7 x 10~^^ecm, using the process e~^e^ 
T+T'j; and by ARGUS [11], which set limits |i?e«)| < 4.6 x lO^^'^ecm and 
\Im{dr)\ < 1.8 X 10~"'^^ecm based on a study of e^e^ t^t^ production. 

In this paper we present the first direct measurement of the r lepton's electric 
dipole moment with a sensitivity in the 10~^''ecm range. We improve on the 
sensitivity of ARGUS by an order of magnitude, by reducing the systematic 
uncertainty in the extraction of rf,-, and by analyzing a much larger sample of 
events. 

We search for CP violating effects at the 7rr vertex in the process e"'"e^ 
7* T^T^ using triple momentum and spin correlation observables. The CP 
violating effective Lagrangian can be expressed as 

Ccp = -idris)fa^''^,Td^A,, (1) 

where the electric dipole form factor dr depends in general on s, the squared 
energy of the r-pair system. In common with other authors, we ignore this 
possible s-dependence, assuming dr{s) = dr, which is constant. {dr{0) cor- 
responds to the electric dipole moment of the r, and we shall use this term 
hereafter.) The squared spin density matrix (A^prod) ^r r-pair production 
e+(p)e-(-p) r+(fc, 5+)r-(-fc, 5_) is given by [12] 

>^prod = -^Im + Mdr)Ml, + Im{dr)M]^ + \dr\''M%, (2) 

MIm = tilkl + + \k\k.pf - S+-S.\k\\l - {k-pf) 

+2(fc-5+)(fc-S_)(|fc|2 + (fco - mrf{k-pf) + 2kl{p-S+){p-S-) 

^ A very strict constraint \dr\ < 2.2 x lO^^^ecm has also been derived from the 
experimental limits on ^ 67 decay [7]. We note, however, that this result as- 
sumes a particular ansatz for the lepton mixing matrix. In a recent preprint [8], 
other authors have argued that the same constraint may be derived under weaker 
assumptions. 
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-2A;o(A;o - mr){k-p){{k-S+){p-S.) + {k-S^){p-S+))], (3) 
M%e = ^Uk\[-{mr + {ko-mr){k-pf){S+xS.)-k + ko{k-p){S+xS.)-p],{^ 

Mjm = ^Uk\[-imr + (fco - mr)ik-pf){S+-S-)-k + koik-p)iS+-S^)-p], (5) 

MI2 = 4e2|fcp-(l - {k-pf)il - (6) 

where ko is the energy of the r, m,- is the r mass, p is the momentum vector 
of e"^, k is the momentum vector of r"^ in the center-of-mass frame, S± are 
the spin vectors for r^, and the hat denotes a unit momentum. We disregard 
the higher order terms proportional to \d'^\, which is vahd since dr is small. 
M 

SM corresponds to the SM term. A^^.^ and Aij^ are the interference terms 
(related to the real and imaginary parts of d^, respectively) between the SM 
and CPV amphtudes. Mj^^ is CP odd and T odd, while Mj^ is CP odd, but 
T even. In the above equations, e'^ and e~ are assumed to be unpolarized and 
massless particles. 

We adapt the so-called optimal observable method [13], which maximizes the 
sensitivity to dr- Here the optimal observables are defined as 



The mean value of the observable Ors is given by 

{One) oc J OneMl^^dcP = J + R<dr) j (8) 

where the integration is over the phase space (0) spanned by the relevant 
kinematic variables. The cross-term containing the integral of the product of 
s-^d Al/m drops out because of their different symmetry properties. The 
expression for the imaginary part is similar. The means of the observables 
{One) and {Oim) are therefore linear functions of dr, 

(One) = ttRe • Re{dr) + hue, {Oim) = aim ■ Im{dr) + him- (9) 



Eight different final states in the decay of r-pairs, {ei'u){fiuu). (ei/p) (tti^) , 
{fii'u)(nu), (euu)(pu), (fiuu)(pu), (7iu)(pu), (pu)(pi?), and (7rz/)(7rz/), are ana- 
lyzed, where all particles except u and u are positively or negatively charged. 
Because of the undetectable particles, we can not fully reconstruct the quan- 
tities k and S±- Therefore, for each event we calculate possible kinematic 
configurations and obtain the mean value of A1sm> -^Re Mjm by av- 
eraging over the calculated configurations. In the case when both r leptons 
decay hadronically (r — > ttu or pu), the r flight direction is calculated with 
a two-fold ambiguity [14] and we take the average of A^sm; At/je and J^jm 
over the two solutions. In the case when either one or both r leptons decay 
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Icptonically (r — > ev9 or i^vP), a Monte Carlo (MC) treatment is adopted 
to take into account the additional ambiguity in the effective mass of the z/z/ 
system {rriup). For each event we generate 100 MC configurations using a hit- 
and-miss approach by varying m^p, and compute the averaged A^sm) •^%e 
and over successful tries in which the r direction can be constructed 

kinematically. In the calculation, we ignore the effect of undetected photons 
coming from initial state radiation, radiative r decays and bremsstrahlung. 
The resulting bias in d^- is included among the systematic errors; the effect on 
the final result is neghgible. 



2 Data and event selection 

In this analysis, we use 26.8 million r-pairs (29.5 fb^^) accumulated with the 
Belle detector [15] at the KEKB accelerator [16]. KEKB is an asymmetric 
energy e+e~ coUider with a beam crossing angle of 22 mrad. Its center-of- 
mass energy is 10.58 GeV, corresponding to the T(45') resonance, with beam 
energies of 8 and 3.5 GeV for electrons and positrons, respectively. Belle is 
a general-purpose detector with an asymmetric structure along the beam di- 
rection. Among the detector elements, the central drift chamber (CDC) and 
the silicon vertex detector (SVD) are essential to obtain the momentum vec- 
tors of charged particles. The combined information from the silica Aerogel 
Cherenkov counters (ACC), the time-of-flight counters (TOF), the Csl elec- 
tromagnetic calorimeter (ECL), and the ii/K^ detector (KLM) is used for 
particle identification. 

The MC event generators KORALB/TAUOLA [17] are used for r-pair produc- 
tion and decays. The detector simulation is performed with a GEANT-based 
program, GSIM. Actual data and MC generated events are reconstructed by 
the same program written by the Belle collaboration. 

We reconstruct the eight r-pair decay modes mentioned above, with the follow- 
ing conditions. Each charged track is required to have a transverse momentum 
Pt > 0.1 GeV/c. Photon candidates should deposit an energy of > 0.1 GeV 
in the EGL. A signal event is required to have two charged tracks with zero 
net-charge and no photon apart from p"^ — > tt^tt^, 7r° — > 77. 

A track is identified as an electron using a likelihood ratio combining dE/dx in 
the GDC, the ratio of energy deposited in the EGL and momentum measured 
in the CDC, the shower shape in the EGL and the hit pattern from the ACC. 
The identification efficiency is estimated to be 92% with a vr^ fake rate of 
0.3% for the momentum range between 1.0 GeV/c and 3.0 GeV/c in the 
laboratory frame [18]. A muon is identified by its range and hit pattern in 
the KLM detector, with efficiency and fake rate estimated to be 91% and 2%, 
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Table 1 



Yield, purity and background rate obtained for the event selection described in the 
text, where the purity is evaluated by MC simulation and its error comes from MC 
statistics. 
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respectively, for momenta in the laboratory frame greater than 1.2 GeV/c. A 
track is considered to be a pion if it is identified as a hadron by the KLM 
information and not identified as an electron: the efficiency of this selection 
is estimated to be 81%, and the purity for the selected samples is estimated 
by MC to be 89%. A p"^ is reconstructed from a charged track and a 7r° 
where the track should be neither an electron nor a muon, and for 7r° — > 77, 
the reconstructed 7r° should have an invariant mass between 110 and 150 
MeV/c^ and a momentum in the laboratory frame larger than 0.2 GeV/c. In 
order to suppress background and to improve the performance of the particle 
identification cuts, we restrict the analysis to lepton candidates within the 
barrel region, —0.60 < cos^ < 0.83, and to pion candidates for r — >■ ttz/ 
within the KLM barrel region, —0.50 < cos^^ < 0.62, where 9 is the polar 
angle opposite to the e+ beam direction in the laboratory frame. For the same 
reason, the laboratory frame particle momentum is required to be greater than 
0.5 GeV/c for an electron, 1.2 GeV/c for both a muon and pion, and 1.0 GeV/c 
for the p"^. 

The dominant backgrounds are due to two-photon as well as Bhabha and iifi 
processes. In order to remove two-photon events, we require the missing mo- 
mentum not to be directed towards the beam-pipe region (imposing a selection 
—0.950 < cos^ < 0.985), and to reject the latter processes we require that the 
sum of the charged track momenta be less than 9 GeV/c in the center-of-mass 
frame. Additional selections are imposed particularly on the eix mode where 
a large number of Bhabha events could contribute through misidentification. 
For the ctt mode, we remove events which satisfy the following criteria: the 
opening angle of the two tracks in the plane perpendicular to the beam axis 
is greater than 175°, and their momentum sum is greater than 6 GeV/c in 
the r-pair rest frame. Finally, we remove events in which the r flight direction 
cannot be kinematically reconstructed, which mostly arise from r-pairs having 
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Fig. 1. Momentum distributions of (a) e , (b) /x , (c) tt , and (d) p in the lab- 
oratory frame. The points with error bars are the data and the histogram is the 
MC expectation. The latter is scaled to the total number of entries. The hatched 
histogram is the background distribution evaluated by MC. 

hard initial-state radiation and misidentified r-pair backgrounds. 

The yield of events passing this selection is given in Table 1 for each of the 
eight selected modes. The mean energy of the r-pair system in the obtained 
sample is ^/s = 10.38 GeV; this sets the scale at which dr{s) is measured. 
Because of events with soft radiated photons, the energy scale is slightly lower 
than the beam energy. The dominant background sources are also listed in 
the table. Hadronic r decays with two or more 7r°'s make a contribution of a 
few percent. For the modes including tt^, kaons and misidentified muons from 
other T decays produce a feed-across background; for example, we estimate 
that 5.3% of the /ivr sample consists of fiK decays. The other backgrounds 
are estimated by MC to be a few percent from two-photon processes, and less 
than 1% from Bhabha, /i/j,, and multihadronic processes. 

Figures 1 and 2 show the resulting momentum and cos 9 distributions, re- 
spectively, for charged particles in the laboratory frame. Very good agreement 
with MC is found, except for low-momentum electrons (Fig. 1(a)) and pions 
(Fig. 1(c)). The dip in the cos^ distribution of the muon (Fig. 2(b)) is due 
to an efficiency drop at the region of overlap between the barrel and endcap 
KLM elements. 

The resulting Ojie and Ojm distributions are shown in Fig. 3 along with those 
obtained from MC simulation with dr — 0. Good agreement is found between 
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Fig. 2. The cos^ distributions of (a) e^, (b) jj,^, (c) ir^, and (d) in the laboratory 
frame. The meanings of the points and histograms are the same as in Fig. 1. 

the experimental data and the MC samples. Distributions of the ratio of the 
data to MC (not shown here) are flat and near 1.0. 



3 Extraction of dr 



In order to extract the value from the observable using Eq. (9), we have to 
know the coefficient a and the offset b. In the ARGUS analysis [11], which also 
used the optimal observable method, the first term of Eq. (8) was assumed to 
be negligible because it vanishes when the integration takes place over the full 
phase space: / M^^dcf) = 0. The value of d^ was obtained as the ratio of the 
observablc's mean to the second term, Re{dr) = {ORe)/{0\^). However, the 
detector acceptance 77 affects the means of the observables according to 

{Ou.) oc j ri{4>)Ou.Ml,^d4^. (10) 

Similar expressions obtain for the imaginary part. This means that the first 
term of Eq. (8) is not necessarily zero and the coefficient may differ from {p\^) 
when the detector acceptance is taken into account. In the ARGUS study, the 
acceptance effect produced the largest systematic uncertainty, of the order of 
10~^^ecm. 

In order to reduce this systematic effect, we extract both parameters a and 
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Fig. 3. Distributions of the optimal observables O^e and Ojm for each mode. The 
upper figure of each mode is for Ors and lower figure is for Ojm- The closed circles 
are the experimental data and the histogram is the MC expectation with dr = 0, 
normalized to the number of entries. The hatched histogram is the background 
distribution evaluated by MC. 



b from the correlation between {ORe){{Oim)) and Re{dr){Im{dr)) obtained 
by a full MC including acceptance effects. An example of the correlation be- 
tween {ORe){{Oim)) and Re{dr){I'm{dr)) is shown in Fig. 4. Each point is 
obtained from MC with detector simulation and event selection. By fitting 
the correlation plot with Eq. (9), the parameters a and b are obtained. The 
feed-across background from other r decays shows some dependence on d^, 
because the spin direction is correlated with the momenta of the final state 
particles. Therefore, the effects of the feed-across background on the coef- 
ficient a and offset b are corrected using the parameters obtained from the 
background. The resulting coefficients and offsets are shown in Fig. 5. 
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Fig. 4. Re[d-j-){Im{dr)) dependence of the mean of the observable {ORe){{Oim)) 
for the Tip mode. The closed circles show the dependence for Re^d-j-) and the open 
circles show the dependence for /m(dT-). The lines show the fitted linear functions. 
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Fig. 5. Sensitivity a (a) and offset b (b) for each mode from MC. The closed cir- 
cles show the parameters for Re{dT-) and the open circles show the parameters for 
Im{dr)- The errors are due to MC statistics. 

As can be seen from the values of the coefficient a, the up and pp modes 
have the highest sensitivities for o?^, while the tttt mode has a somewhat lower 
sensitivity. For the real part, this is an effect of averaging over the two solutions 
for the r direction: in the r — > tti/ case, this causes spin correlation information 
to be lost, whereas for t ^ pv the angular distribution of the p titi^ decay 
provides information on the r spin which survives the averaging procedure. 
For the imaginary part, the lower sensitivity is due to the tighter cos 6' cut 
applied to pions in the r — > tti/ channel, compared to r ^ pv] there is also 
a small effect from the tighter momentum cut. The remaining modes, which 
include leptons, have low sensitivities because information about the r spin, 
and also its direction, is lost due to the additional neutrino (s). Non-zero offsets 
bim are seen for the imaginary part, due to the forward/backward asymmetry 
in the acceptance of the detector. 
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4 Systematic uncertainties 



Although in general the MC simulation reproduces the observed kinematic 
distributions well, the small disagreement that is evident in Figures 1 and 2 

dominates the systematic uncertainty. The effect on d^. is studied by reweight- 
ing the MC distributions by the ratio of data and MC. The second significant 
uncertainty originates from possible charge asymmetry in the detection effi- 
ciency. The ratio of yields, N [a^ I3~) / N {a~ (3^) , for data and MC agrees within 
1%, where a and (3 are the relevant charged particles from the r decays. The re- 
sulting systematic uncertainty is evaluated by varying the detection efficiency 
by ±1%. The effect is of the same size as the statistical error for Im{dT-), 
while it is negligible for Re{dT^). The backgrounds lead to additional system- 
atic uncertainties in dr because the parameters a and h are corrected for the 
background distributions: this effect is assessed by varying the assumed back- 
ground rate. The effects of photon energy resolution, and of possible biases 
in reconstructed momentum (for charged tracks and photons), are checked by 
applying scaling factors based on a comparison of data and MC distributions. 
In order to examine a possible asymmetry arising from the alignment of the 
tracking devices, we measure the differences in polar angles, A9 — ^Jm ~ ^cm) 
and azimuthal angles, A(j) = ^^j^ — 0cm > the two tracks in e+e" l^'^IJ'~ 
events, and find a small deviation from back-to-back topological alignment in 
each direction: A6 = 1.48 mrad and A(p = 0.36 mrad. Applying an artificial 
angular deviation of this magnitude to one of the charged tracks, we find the 
change in the observables to be negligible compared to the other errors. 

To estimate the effect of ignoring radiative processes in the calculation of the 
observables, we consider the case dr = and compare our standard calculation 
to one which includes initial state radiation, taking the shift in d,- as a measure 
of the systematic error. Since this shift occurs in analysis of both data and 
MC events, it is already taken into account in the analysis (up to effects of 
detector and/or background mismodelling), so the estimate is conservative. 
For all decay modes apart from tttt, the shift is negligibly small. 

The various sources of systematic error are hsted in Table 2. 



5 Result 

The values of dj- extracted using Eq. (9) are hsted in Table 3 along with the 
corresponding statistical and systematic errors, and plotted in Fig. 6. The 
results are consistent with d-^ — within the errors. 

Finally, we obtain mean values for Re{dr) and /m(d^) over the eight different 
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Table 2 

Systematic errors for Re{dr) and Im{dr) in units of 10 ^^ecm. 
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0^ 
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10 
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Al 


Charge asymmetry 


0.13 


0.44 


0.43 


0.02 


0.09 


0.15 






Background variation 


0.08 


0.07 


0.02 


0.01 


0.03 


0.02 


0.03 


0.06 


Momentum reconstruction 


0.03 


0.03 


0.06 


0.00 


0.02 


0.02 


0.04 


0.04 


Detector alignment 


0.01 


0.03 


0.02 


0.01 


0.01 


0.02 


0.01 


0.05 


Radiative effects 


0.05 


0.02 


0.01 


0.01 


0.01 


0.02 


0.01 


0.02 


Total 


0.46 


0.45 


0.44 


0.13 


0.10 


0.16 


0.11 


0.42 



Tabic 3 

Results for the electric dipole moment. The first error is statistical and the second 
is systematic^ 



Mode 


Re{dr) (lO-^^ecm) 


Im{dr) (lO-i^ecm) 


eji 


2.25 ± 1.26 ± 0.93 


-0.41 ± 0.22 ± 0.46 




0.43 ± 0.64 ± 0.60 


-0.22 ±0.19 ±0.45 


/XTT 


-0.41 ±0.87 ±0.74 


0.15 ±0.19 ±0.44 


ep 


0.00 ±0.36 ±0.14 


-0.01 ±0.14 ±0.13 


HP 


0.04 ±0.42 ±0.18 


-0.02 ±0.14 ±0.10 


np 


0.34 ±0.25 ±0.22 


-0.22 ±0.13 ±0.16 


PP 


-0.08 ±0.25 ±0.17 


-0.12 ±0.14 ±0.11 


TT 77 


0.12 ± 1.17 ± 0.18 


0.2-1 ± 0.31 ± 0.-12 


Mean value 


0.115 ±0.170 


-0.083 ± 0.086 
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Fig. 6. Re{dr) and Im{dr) for each mode. Both statistical and systematic errors are 
included. The small ticks on the error bars show the statistical errors. 

T+T~ modes weighted by quadratically summed statistical and systematic 
errors, 

Re{d^) = (1.15 ± 1.70) X 10-^^ecm, (11) 
Im{dr) = (-0.83 ±0.86) x 10-^^ecm, (12) 

with corresponding 95% confidence limits 

- 2.2 < Re{dr) < 4.5 (10~^^ecm), (13) 
-2.5 < 7m(d^) < 0.8 (lO'^^ecm). (14) 



This investigation has improved the sensitivity to the r lepton's electric dipole 
moment by an order of magnitude over previous measurements. 



Acknowledgements 



We would like to thank Professors K. Hagiwara, O. Nachtmann, and Z. W§,s for 
their constructive advice and many helpful discussions. We also wish to thank 
the KEKB accelerator group for the excellent operation of the KEKB acceler- 
ator. We gratefully acknowledge the support from the Ministry of Education, 
Culture, Sports, Science, and Technology of Japan, Grant-in-Aid for JSPS 
Fellows 01655 2001, and the Japan Society for the Promotion of Science; the 
Australian Research Council and the Australian Department of Industry, Sci- 
ence and Resources; the National Science Foundation of China under contract 
No. 10175071; the Department of Science and Technology of India; the BK21 
program of the Ministry of Education of Korea and the CHEP SRC program 
of the Korea Science and Engineering Foundation; the Polish State Commit- 
tee for Scientific Research under contract No. 2P03B 17017; the Ministry of 
Science and Technology of the Russian Federation; the Ministry of Education, 
Science and Sport of the Republic of Slovenia; the National Science Council 
and the Ministry of Education of Taiwan; and the U.S. Department of Energy. 



14 



References 



[I] K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 87 (2001) 091802; Phys. 
Rev. D 66 (2002) 032007. 

[2] B. Aubert et al. (BaBar Collaboration), Phys. Rev. Lett. 87 (2001) 091801; 
Phys. Rev. D 66 (2002) 032003. 

[3] T. Huang, W. Lu, and Z. Tao, Phys. Rev. D 55 (1997) 1643. 

[4] W. Bernreuther, A. Brandenburg, and P. Overmann, Phys. Lett. B 391 (1997) 
413; ibid., B 412 (1997) 425. 

[5] K. Akama, T. Hattori, and K. Katsuura, Phys. Rev. Lett. 88 (2002) 201601. 

[6] R. Escribano and E. Masso, Phys. Lett. B 301 (1993) 419; Nucl. Phys. B 429 
(1994) 19; Phys. Lett. B 395 (1997) 369. 

[7] T. Huang, Z. H. Lin, and X. Zhang, Phys. Lett. B 450 (1999) 257. 

[8] J. Q. Zhang, X. C. Song, W. J. Huo, and T. F. Feng, |hep-ph/020530£ . 

[9] M. Acciarri et al. (L3 Collaboration), Phys. Lett. B 434 (1998) 169. 

[10] K. Ackerstaff et al. (OPAL Collaboration), Phys. Lett. B 431 (1998) 188. 

[II] H. Albrecht et al. (ARGUS Collaboration), Phys. Lett. B 485 (2000) 37. 

[12] W. Bernreuther, O. Nachtmann, and P. Overmann, Phys. Rev. D 48 (1993) 78. 

[13] D. Atwood and A. Soni, Phys. Rev. D 45 (1992) 2405. 

[14] K. Ackerstaff et al. (OPAL Collaboration), Z. Phys. C 74 (1997) 403. 

[15] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Meth. A 479 (2002) 
117. 

[16] E. Kikutani ed., KEK Preprint 2001-157 (2001), to appear in Nucl. Instr. and 
Meth. A. 

[17] KORALB(v.2.4)/TAUOLA(v.2.6): S. Jadach and Z. W^s, Comp. Phys. 
Commun. 85 (1995) 453; ibid., 64 (1991) 267; ibid., 36 (1985) 191; S. Jadach, 
Z. Wcis, R. Decker, and J.H. Kiihn, Comp. Phys. Commun. 64 (1991) 275; ibid., 
70 (1992) 69; ibid., 76 (1993) 361. 

[18] K. Hanagaki et al., Nucl. Instr. and Meth. A 485 (2002) 490. 



15 



